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We construct density dependent Class /// charge symmetry violating (CSV) potential due to 
mixing of p-co meson with ofT-shell corrections. Here in addition to the usual vacuum contribution, 
the matter induced mixing of p-uj is also included. It is observed that the contribution of density 
dependent CSV potential is comparable to that of the vacuum contribution. 
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I. INTRODUCTION 

The exploration of symmetries and their breaking have 
always been an active and interesting area of research in 
nuclear physics. One of the well known examples, that 
can be cited here, is the nuclear /? decay which violates 
parity that led to the discovery of the weak interaction. 
Our present concern, however, is the strong interaction 
where, in particular, we focus attention to the charge 
symmetry violation (CSV) in nucleon-nucleon (NN) in- 
teraction. 

Charge symmetry implies invariance of the NN inter- 
action under rotation in isospin space, which in nature, 
is violated. The CSV, at the fundamental level is caused 
by the finite mass difference between up (u) and down 
(d) quarks [J, H, d, |3j H, 0|- As a consequence, at the 
hadronic level, charge symmetry (CS) is violated due to 
non-degenerate mass of hadrons of the same isospin mul- 
tiplet. The general goal of the research in this area is 
to find small but observable effects of CSV which might 
provide significant insight into the strong interaction dy- 
namics. 

There are several experimental data which indicate 
CSV in NN interaction. For instance, the difference 
between pp and nn scattering lengths at ^Sq state is 
non-zero 0, 0i Si- Other convincing evidence of CSV 
comes from the binding energy difference of mirror nu- 
clei which is known as Okamoto-Nolen-Schifer (ONS) 
anomaly i, El- 

The modern manifestation of 
CSV includes difference of neutron-proton form factors, 
hadronic correction to g — 2 the observation of the 
decay of ^''(3686) ^ {J/^)Tr° etc [13]. 

In nuclear physics, one constructs CSV potential to 
see its consequences on various observables. The con- 
struction of CSV potential involves evaluation of the 
NN scattering diagrams with intermediate states that 
include mixing of various isospin states like p-u or tt- 
77 mesons. The former is found to be most dominant 
11, [H, El, [H, m, [13 which we consider here. 

Most of the calculations performed initially to con- 
struct CSV potential considered the on-shell [l^ or con- 
stant p-uj mixing amplitude [l7j , which are claimed to be 
successful in explaining various CSV observabl es [itI . [isj . 
This success has been called into question [l^, |2C| on the 
ground of the use of on-shell mixing amplitude for the 



construction of CSV potential. First in [20| and then 
in [m, [H [H, [H [H], it is shown that the p-uj mixing 
has strong momentum dependence which even changes 
its sign as one moves away from the p (or w) pole to the 
space-like region which is relevant for the construction of 
the CSV potential. Therefore inclusion of off-shell correc- 
tions are necessary for the calculation of CSV potential. 
We here deal with such mixing amplitude induced by the 
N-N loop incorporating off-shell corrections. 

In vacuum, the charge symmetry is broken explicitly 
due to the non-degenerate nucleon masses. In matter, 
there can be another source of symmetry breaking if the 
ground state contains unequal number of neutrons (n) 
and protons (p) giving rise to ground state induced mix- 
ing of various charged states like p-uj meson even in the 
limit M„ = Mp. This additional source of symmetry 
breaking for the construction of CSV potential has, to 
the best of our knowledge, not been considered before. 

The possibility of such matter induced mixing was 
first studied in [26| and was subsequently studied in 
[27I . [H, [2^, [3^ . For the case of Tr-ry meson also such asym- 
metry driven mixing is studied in [3l[ . But none of these 
deal with the construction of two-body potential and the 
calculations are mostly confined to the time-like region 
where the main motivation is to investigate the role of 
such matter induced mixing on the dilepton spectrum 
observed in heavy ion collisions, pion form factor, meson 
dispersion relations etc. [13, [lol- In Ref . [3^ . attempt 
has been made to calculate the density dependent CSV 
potential where only the effect of the scalar mean field 
on the nucleon mass is considered excluding the possibil- 
ity of matter driven mixing. All existing matter induced 
mixing calculations, however, suggest that, at least in the 
p-uj sector, the inclusion of such a matter induced mix- 
ing amplitude into the two body NN interaction poten- 
tial can significantly change the results both qualitatively 
and quantitatively. It is also to be noted that such mix- 
ing amplitudes, in asymmetric nuclear matter (ANM), 
have non-zero contribution even if the quark or nucleon 
masses are taken to be equal [H, [13, [2!, [H [s^]. We 
consider both of these mechanisms to construct the CSV 
potential. 

Physically, in dense system, intermediate mesons 
might be absorbed and re-emitted from the Fermi 
spheres. In symmetric nuclear matter (SNM) the emis- 
sion and absorption involving different isospin states like 



2 



' vac 



rNN 

m.ed 



FIG. 1: Feynman Diagrams that contribute to the construc- 
tion of CSV NN potential in matter. Solid lines represent 
nucleons and wavy lines stand for vector mesons. 



p and LU cancel when the contributions of both the pro- 
ton and neutron Fermi spheres arc added provided the 
nucleon masses are taken to be equal. In ANM, on the 
other hand, the unbalanced contributions coming from 
the scattering of neutron and proton Fermi spheres, lead 
to the mixing which depends both on the density (ps) 
and the asymmetry parameter [a = {pn — Pp)/ ps]- In- 
clusion of this process is depicted by the second diagram 
in FigU] represented by which is non-zero even in 

symmetric nuclear matter if explicit mass differences of 
nucleons are retained. In the first diagram, involves 
NN loop denoted by the circle. The other important el- 
ement which we include here is the contribution coming 
from the external legs. This is another source of explicit 
symmetry violation which significantly modify the CSV 
potential in vacuum as has been shown only recently by 
the present authors [33| . 

This paper is organized as follows. In Sec. II we present 
the formalism where the three momentum dependent p^- 
uj mixing amplitude is calculated to construct the CSV 
potential in matter. The numerical results are discussed 
in Sec. III. Finally, we summarize in Sec. IV. 



II. FORMALISM 

We start with the following effective Lagrangians to 
describe ujNN and pNN interactions: 



jNN 



c 



pNN 



7m 



(la) 

f.#^vlr, (lb) 



where Cp — fp/gp is the ratio of vector to tensor 
coupling, M is the average nucleon mass and f is the 
isospin operator. and represent the nucleon and 
meson fields, respectively, and g's stand for the meson- 
nucleon coupling constants. The tensor coupling of to is 
not included in the present calculation as it is negligible 
compared to the vector coupling. 

The matrix element, which is required for the construc- 
tion of CSV NN potential is obtained from the relevant 
Feynman diagram [33i |: 



= [tx^(p3)r^(g)«^(Pi)] A^„(g)n^5(g2) 

X A^,(g) [uN{pi)T:,{~q)uN{p2)]. (2) 



In the limit qo — > 0, Eq.Q gives the momentum space 
CSV NN potential, V^gvi^)- Here r£;(g) = g^7^, 
Tp{q) = gp 7*^ — ^ijif^^'^lx denote the vertex factors, 

Wat is the Dirac spinor and A^^(q), (i = p,uj) is the me- 
son propagator, pj and q are the four momenta of nucleon 
and meson, respectively. 

In the present calculation, p-u mixing amplitude (i.e. 
polarization tensor) n^^((7^) is generated by the differ- 
ence between proton and neutron loop contributions: 

If;;:iq') ^Il';:^P\q') ~U';:^-\q'). (3) 
Explicitly, the polarization tensor is given by 

fflJ^:^^'(9')-/(0iTr[r^;(g)G^(fc)r;;(-g)GA.(fc+g)], 

(4) 

where k = (fco,k) denotes the four momentum of the 
nucleon in the loops (see Fig[T|) and Gat is the in- medium 
nucleon propagator consisting of free (G^) and density 
dependent (G^) parts [33 |. 



G^(fc) 



(5a) 



G^(fc) 



N 



-(^ + MN)S{ko - EN)9{kN - |k|) (5b) 



The subscript N stands for nucleon index (i.e. N ~ p ov 
n), denotes the Fermi momentum of nucleon, nucleon 
energy En = ^/A/^T+^w nucleon mass is denoted 
by Mat. 9{kN ~ |k|) is the Fermi distribution function at 
zero temperature. 

The origin of G^(fc), in addition to the free propagator 
resides in the fact that here one deals with vacuum con- 
taining real particles which when acted upon the annihi- 
lation operator does not vanish (see Appendix A for de- 
tails). The appearance of the delta function in Eq. (j5bp in- 
dicates the nucleons are on-shell while d{kN — \k\) ensures 
that propagating nucleons have momentum less than k^. 

Likewise, the polarization tensor of Eq.(|4]) also con- 
tains a vacuum [Ilvac^\q'^)] and a density dependent 
Pm'ld^'' (9^)] parts as shown in Fig[TJ It is to be noted 
that the density dependent part given by the combina- 
tion of GfjG^ + G^G^ corresponds to scattering that we 
have discussed already, whereas the term proportional to 
G^G^ vanishes for low energy excitation [33. The vac- 
uum part, viz. [Ilvac^\q'^)] on the other hand involves 
G^G^ which gives rise to usual CSV part of the poten- 
tial due to the splitting of the neutron and proton mass. 
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It might be worthwhile to mention here that Eq. (|5b[) 
can induce charge symmetry breaking in asymmetric nu- 
clear matter due to the appearance of the Fermi distri- 
bution function in the propagator itself which can dis- 
tinguish between neutron and proton, even if their mass 
are taken to be degenerate. Evidently, this is an exclu- 
sive medium driven effect where, as mentioned already in 
the introduction, the charge symmetry is broken by the 
ground state. The total charge symmetry breaking would 
involve both the contributions where, it is clear that even 
for a = 0, the medium dependent term can contribute if 
the non-dcgcncratc nuclcon masses are considered. 

Note that the polarization tensor ^p^{q^) can be ex- 
pressed as the sum of longitudinal component ^^^^(g^)] 
and transverse component [nj^((7^)] which will be useful 
to simplify the matrix element given in Eq.Q. 



that both n^^^(g2) and n^''erf((?^) obey the current con- 
servation as Qf.Q'"' = q^Qf"" = and Qf^K'"' = qt^Kf"" = 
0. The superscripts vv and tv in Eqs.lH]) - (fTO|) indi- 
cate the vector- vector and tensor- vector interactions, re- 
spectively. The dimensional counting shows that vac- 
uum part of the polarization tensor [both Ilvac^\q'^) and 
n'a'c^'' (g^)] is ultraviolet divergent and dimensional reg- 
ularization [s^, [H, [s^ is used to isolate the divergent 
parts. Since the mixing amplitude is generated by the 
difference between the proton and neutron loop contribu- 
tions, the divergent parts cancel out yielding the vacuum 
amplitude finite. 



n^/Jq') ==Tl^,Jq')A^'' +TlUq')B^''' 



(6) 



where A'"' and i?'"' are the longitudinal and transverse 
projection operators [H]. We define 11^^ = -11^° + 11^3 

and nj^ = nii = n22. 

We, in the present calculation, use the average of lon- 
gitudinal and transverse components of the polarization 
tensor instead of 11^^ and 11 J^. The average mixing am- 
plitude is denoted by 



^q^Q^" / dx 



X In 



x)q 



M'^ — x(l — 



(11) 



fl(g2) 



iq') 
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3 

= nuaciq^)+fim.ediq'^)- 



(7) 



In the last line of Eq.([7]), ^vac{q^) and Iimed{q^) denote 
the average mixing amplitudes of vacuum and density 
dependent parts, respectively. 

To obtain tVvac{q^) and tVmed{q^) one would calculate 
the total polarization tensor given in Eq.(l4]). After eval- 
uating the trace of Eq.Q, we find the following vacuum 
and density dependent parts of the polarization tensor. 



Eq. pip shows the four- momentum dependent vacuum 
polarization tensor. From the above equation one can 
calculate longitudinal (n^^jc) ^'^d transverse (H^^^) com- 
ponents of the vacuum mixing amplitude and in the limit 
go ^ 0, n^^^(q^) = n^^^(q^). Therefore, the average 
vacuum mixing amplitude is 



^Zn^l') = Q^''\Ktf\q') + KL^\q')\. (8) 



and 



n 



Hv{vv{N)} ( 2 



{2'Kf2EN 
q^Rf" - {q ■ kfQf"' 



-3777^e(fc^ - |k|) 



g4 - A{q ■ kf 

nrir^^»(g^) = 4g,g.C, / - |k|) 



(9) 



{2ttY2En 



q^Q 



g4 - A{q ■ fc)2 



(10) 



where Q^''' = {-gf" + q''q''/q'^) and Kf"" 
(^kt" - {q ■ k)^'j (^k" - {q ■ k)^y It is to be mentioned 



n™c(q') = -[ni,(q2) + 2n^,,(q2)] 



127r2 



(2 + 3Cp) In 



Mp 



(12) 



Eq. (|T2)l represents the three momentum dependent 
vacuum mixing amplitude. This mixing amplitude van- 
ishes for Af„ — Mp and then no CSV potential in vacuum 
will exist. 

To calculate density dependent mixing amplitude from 
Eq.® and (HH]) we consider « Mjy. In the hmit 
qo — *■ 0, one finds following expressions: 
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(13) 
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q — 2fcAr 
q+ 2kN 



(14) 



Note that the terms within the first curly braces of 
both Eg. p^ and Eg. p^ arise from the vector- vector 
interaction while the terms within the second curly 
braces arise from tensor-vector interaction of the den- 
sity dependent polarization tensor. The 33 component 
of the density dependent polarization tensor vanishes i.e. 
n^l^^(q^)=0. Now 



where 



n,„ed(q^) 

nLd(q') 
nLd(q') 



3 [n^ed(q') + 2nLd(q') 



^^med vq ) ^^med Vq ) 
^'^med vq / ^'^med ) 



(15) 

(16) 
(17) 



With the suitable expansion of Eqs. (fl3|) and ([T4|) in 



terms of |q|/fcp(„) and keeping 0(q /fcp(„) 
nw(q2) ^ A'-^'q2, 



where 
A' = 

A' 
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terms wc get 
(18) 
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(20) 



Clearly, limed (q^) is also three momentum dependent 
and if Af„ = Mp it vanishes in SNM but is non- vanishing 
in ANM. In the present calculation nucleon masses are 
considered non-degenerate and the asymmetry parameter 



To construct CSV potential we take non-relativistic 
(NR) limit of the Dirac spinors in which case we obtain 
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(21) 



where a is the spin of nucleon. P denotes the average 
three momentum of the interacting nucleon pair and q 
stands for the three momentum of the meson. 

The explicit expression of full CSV potential in mo- 
mentum space can be obtained from Eq.(13) of ref.fssj 
by replacing the mixing amplitude Ilpi^{q) with Il(q^). 

gpgu n(q2) 



(q2 +m2)(q2 4-to2- 



T. 



3P2 



2Ml Mil AMI 
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4A4 
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^ (cTi ■ q)(cr2 • q) 
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^ 2A/ \ V2A/ 2A/^ ' 

1 , w AA/(1,2) 
— (.,.q)(a2.q)^^ 



^(o-i -CT2) • (qx P) 



(22) 



Eq. 



presents the full CSV NN potential in mo- 



mentum space in matter. Here Tg = r3(l) ± r3(2) and 
S = ^(ui -1-172) is the total spin of the interacting nucleon 
pair. We define M = (A/„-|-A/p)/2, AAf ^ (A-/„-A.fp)/2 
and AAf(l,2) = -AA/(2, 1) = AM. It is be men- 
tioned that the spin dependent parts of the potential 
appear because of the contribution of the external nu- 
cleon legs shown in Fig[TJ On the other hand, 3P2/2Af^ 
and — q2/8Af^ arise due to expansion of the relativistic 
energy Epf of the spinors. 

In matter, Vf^^{q), consists of two parts, one con- 
tains the vacuum mixing amplitude and other contains 
the density dependent mixing amplitude. The former 
is denoted by ^-^^^(q) and later by V,^J^{c[). Thus, 

^c''5^(q) = Kr(q) + v;S(q)- 

From Eq. (P^ we extract the following term which, in 
coordinate space gives rise to the J-function potential. 



^yNN 

O^CSV 



'.pguj (A + 
1 + 2C, 
8A4 



A') 



1 + Cp 
AMj, 



(di • (T2) 



.(23) 



To avoid the appearance of (5-function potential in 
coordinate space, one should introduce form factors 
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[Fi(q^)], for which the meson- nucleon coupHng constants 



become momentum dependent i.e. gi(q^) 
Here we use the following form factor. 



A? 



A? 



(24) 



where A,; is the cut-off parameter governing the range of 
the suppression and m,; denotes the mass of exchanged 
meson. 



The full CSV potential presented in Eq. (|22|) contains 
both Class /// and Class IV potentials, and both break 
the charge symmetry in NN interactions. The terms 
within the first and the second curly braces represent 
Class (///) and Class IV potentials, respectively. Class 
(///) potential differentiates between nn and pp systems 
while Class (IV) NN potential exists in the np system 
only. We, in this paper, restrict ourselves to Class (///) 
potential only. 

The coordinate space potential can be easily obtained 
by Fourier transformation of VQgy{q) i.e. 
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^ ( mlYo{xp)-mlYQ{x^) 
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Ml 
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mi 



- mi 




(27) 



where Xi = miV. The explicit expressions of Vvv{xi) and 
Vty{x{) are given in Ref.lH. In Eqs.lUS]) and (27]), the 
M^"^ independent terms represent central parts without 
contributions of external nucleon legs. 
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We drop the term ZP^ /2Ml from Eq.(1221) while taking 
the Fourier transform as it is not important in the present 
context. However, this term becomes important to fit ^Sq 
and phase shifts simultaneously. 

Now the CSV potential in coordinate space without 
^^csv reduces to 
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(26) 



and 



The potentials presented in Eas. ((26l) and (|27p do not 
include the form factors so that these potentials di- 
verge near the core. The problem of divergence near 
the core can be removed by incorporating form factors 
as discussed before. With the inclusion of form factors, 
Vjl^{r) and V^J^{r) take the following form: 



V, 
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[ m2 - m2 
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Ml y ml - ml J 
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In Eqs.dlHl) and dM]), X, = A^r and 



(29) 



parameter a= 1/3. The dashed and dotted curves show 
AT^ac and ISVmed, respectively. The total contribution 
{i.e. IXVvac + lAVmed) is shown by the soHd curve. It is 
observed that the density dependent CSV potential can 
not be neglected while estimating CSV observables such 
as binding energy difference of mirror nuclei. 
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FIG. 2: AV for ^So state (without SVci\/ and form factors) 
at pB = 0.148 fm"^' and a = 1/3. 
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(30) 
(31) 



Note that Eqs. ((28|) and (j29|) contain the contribution 
of SV^gy and the problem of divergence near the core is 
removed. 



III. RESULTS 



FIG. 3: Ay for (with form factors including the Fourier trans- 
form of Wcsv) ^So state at ps = 0.148 Im'^* and a = 1/3. 



Fig [3] displays AV for ^Sq state including the Fourier 



transform of SV^^y and form factors. 



^V,,„.r and AI4. 



are represented by dashed and dotted curves. Solid curve 
shows the sum of these two contributions. Note that in- 
corporation of form factors remove the problem of diver- 
gence near the core. It is observed that the CSV 
potential changes its sign due to inclusion of the Fourier 
transform of <5V(5gy. 



Using Eqs. (p6)) and (|27p we show the difference 
of CSV potentials between nn and systems i.e. 
Ay = V^-^y - Vl\y in Figllfor ^So state at nuclear 
matter density pB = 0.148 fm^"^ with asymmetry 



IV. SUMMARY AND DISCUSSION 

In this work we have constructed the CSV NN po- 
tential in dense matter using asymmetry driven momen- 
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turn dependent p'^-w mixing amplitude within the frame- 
work of one-boson exchange model. Furthermore, the 
correction to the central part of the CSV potential due 
to external nucleon legs are also considered. The closed- 
form analytic expressions both for vacuum and density 
dependent CSV NN potentials in coordinate space are 
presented. 

We have shown that the vacuum mixing amplitude and 
the density dependent mixing amplitude are of similar 
order of magnitude and both contribute with the same 
sign to the CSV potential. The contribution of density 
dependent CSV potential is not negligible in comparison 
to the vacuum CSV potential. 



bks\'^Q) 
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for all k 
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for |k < fc^v 




= n(k)|*o) 





n(k) is either or 1 and this can be accomplished with 
the step function 6{kN - |k|). Using Eqs. (|^3)) -( (X6|) one 
obtains 



APPENDIX A 

The position space Fermion propagator in vacuum is 
given by the vacuum expectation value of the time or- 
dered product of Fermion fields. 



iGN{x-x') = {Q\T[iP{x)'4){x')]\{)) . 



(Al) 



In medium, the vacuum |0} is replaced by the ground 
state ) which contains positive-energy particles with 
same Fermi momentum fcjv and no antiparticles. Thus, 



iGN{x-x') = (*o|^(^)V5(a:')l*o)^(i"i') 

- {■^S{x')^{x)\^^)e{t' ~t). (A2) 



Note that the time-ordered product in Eq. (jA2p in- 
volves negative sign for Fermions. The Fermion field con- 
tains both the positive- and negative-energy solutions. 
The modal expansion for the Fermion fields are. 
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Here and aks arc the creation and annihilation op- 
erators for particles and likewise and bks are the cre- 
ation and annihilation operators for antiparticles. The 
only nonvanishing anticommutation relations are 



{flfe., 4's'} = {bks.biA - 5,s'5\\^ - k'). (A5) 

Since contains only positive-energy particles, we 
have the following relations: 
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and 



(^'o|^(x')^(x)|*o) = 

^[(*0|4's'«fc«l*0) 



d^k 



d^k' 



^{2^f2EkJ ^/WF2Ey 



Uk's'Ukse 



— i{k-x—k'-x ) 



+ ( ^^\hk's'hU^^ ) Vk's'Vkse'^"-''-'''--''^ 



{2TTf2Ek 



(A8) 



Now, 



9{t-t')e 



i\ —ik-{x—x') 



I / —. -. (A9) 

dfco e-*'=-(=^-=^') 



27r fcn - El 



- iifc + «e 



27r ko — Ek — ie 



, (AlO) 



r rlk-„ pik-{x-x) 

e(t'-Oe-(— — . (All) 



J 2n ko-Ek + ie 
From Eqs.dMI, (!A10|) and (jAlip we have 

-W.^e-(-')(^ + M.)-^(^--|^") 



i2n)^2Ek 
d-^k 



ko - Ek- ie 



{2^Y2Ek 



- Mn) 



pik-{x — x ) 

ko - Ek + ie' 



(A12) 



Now changing k —k in the last integral of Eq. (|A12[ ) 
and substituting Eqs.jSI, (EH) and ((XT2|l in Eq. ((X2|) 
we get, 
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From Eqs. (|Il3)) - ((JT5l) . 



■r ( M ■ f ^ 



1 



kf) — EN+ie ko — En —ie kt^ + E^ — ie 



(A13) 



In Ea. (|A13|) E^ has been replaced by E^- The first 
term of Ea. (|A13p represents partiele propagation above 
the Fermi sea and the second term indicates the propa- 
gation of holes inside the Fermi sea. The last term shows 
the propagation of holes in the infinite Dirac sea. Now. 



1 



1 



2E 



N 



ko — EN + if- ko + E^-ie k'^-M'^ + iC 



(A14) 



1 



1 



k^ — En ~ie ko — EN + ie 



= 2iTTS{ko-EN). (A15) 



iGN{x — x') = i 



df^k 
(2^ 



-ik-( 



-^^GN{k) (A16) 



where GN{k) is the in- medium Fermion propagator in 
momentum space. 



_ HMn 



k^-AI^ + ie 



+ {^ + MN)S{ko-EN)e{kN-\k\) 
En 

= GUk) + G^{k). (A17) 
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